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We present a model of leptoquarks (LQs) with a significant partial branching ratio into an extra sector, taken 
to be a viable dark matter candidate, other than the canonical lepton and jets final state. For LQs with mass 
around 500 GeV, the model reproduces the recent excess claimed by the CMS collaboration in jjpv final 
state: the event rate, the distribution in di-lepton invariant mass and the rapidity range are compatible with the 
data. The model is compatible with other collider bounds including LQ searches, as well as bounds from meson 
mixing and decays. Prospects of discovery at Run II of the LHC are discussed. 


Introduction. The CMS collaboration reported a 2.6 ct 
excess compared with Standard Model expectations in 
jj]f}y events, containing two opposite-sign same-flavor 
leptons = {e^, /r^}, at least two jets and missing trans¬ 
verse momentum J/if |1T] with 19.4 fb~^ of integrated luminos¬ 
ity at a center of mass energy of 8 TeV. The excess was found 
in the central region with lepton pseudo-rapidities \rji\ < 1.4, 
after event selection and flavor subtraction cuts, and with di¬ 
lepton invariant mass mu < 80 GeV, as shown in Fig. [T] No 
excess is seen in other regions nor in the trilepton channel. 

The excesses were found in the context of searches for 
edges in mu- The triangular edge is a classic supersymme¬ 
try signal. Interpretations of the CMS excess in the context 
of the so-called golden cascade (x2 -A -A 

have been proposed (x?, X 2 > ^ '^^e lightest neutralino, 

the next-to-lightest neutralino and the slepton, respectively). 
Since direct electroweak production of X 2 has a too small 
cross section to provide a large enough rate whilst evading 
previous collider bounds from LEP, assistance from colored 
particle production is required. The decay chain could start 
with i —t tx 2 12 or q qx^ |f3| . The former interpretation is 
constrained by the fact that the CMS study did not observe a 
large excess in trilepton final states, which should be present 
from the leptonic top quark decay. The CMS study itself opted 
for an explanation in terms of light sbottoms b —t 6 x 2 , and 
X2 -A although no 6-jet requirement was made on 

the final states. Ref. a explored the parameter space in order 
to simultaneously satisfy bounds from 4 charged lepton pro¬ 
duction. Recently, however, several potential SUSY scenarios 
explaining the CMS excess have been shown to be in tension 
with existing experimental data fSl. 

The purpose of this letter is to point out that the CMS excess 
can be explained by a different, non-supersymmetric, class of 
models with leptoquarks (LQs) ii-illl. In a model of LQs 
presented in ifTSl by some of the authors, the LQ had branch¬ 
ings to three possible final states: (i) the usual final state with 
charged lepton and jet; (ii) final state with jets and missing 
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FIG. 1. Di-lepton invariant mass spectrum of CMS £'^£~jjpT 
events: expected non-DY background (black histogram), expected 
DY (red histogram), expected signal for benchmark C (green his¬ 
togram) and expected signal plus background for benchmark C (blue 
histogram).21 signal events are expected for mu > 100 GeV, which 
is compatible with CMS data. 


energy py in the form of dark matter (DM); and (Hi) final 
state with charged lepton, jets, and py in the form of DM. 
By a combination of branchings and kinematics, the model 
was able to account for the mild excess in the eejj and eivjj 
channels observed by a recent CMS search IS El for first gen¬ 
eration LQs in the mass range 550 - 650 GeV. To account 
for the event rate observed by CMS in these first generation 
LQ searches, the branching ratio of the LQ into electrons and 
quarks was taken to be « 15%. A viable DM sector was con¬ 
structed which accounted for the remaining branching ratio 
through channels {ii) and {Hi) mentioned above. 

In this work, we note that the same framework, through 
its connection to DM, interestingly enables us to fit the CMS 
£~^£~jjpy excess. Fig. [^exemplifies how the £'^£~ jjfry ex¬ 
cess can be reproduced by LQs with « 500 GeV masses: the 
model predicts a peak in mu, which fits data roughly as well 
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as the triangular edge searched for by CMS. The model also 
predicts small numbers of events in the forward region and 
for trilepton final states, also compatible with the CMS mea¬ 
surements. The model is compatible with flavor constraints, 
as well as constraints from other LHC searches. We will see 
that the benchmark point that best fits the CMS jjpY 
excess is around ^ 500 GeV where there is no excess in the 
eejj and evjj channels. This necessitates a different choice 
of branchings and benchmark from iia. 

The Model. We consider LQs that do not lead to pro¬ 
ton decay at renormalizable level; either a scalar in the 
(3, 2, 7/6) representation of SU(3) 0 SU(2) ® U(l) that can 
couple to Qe and to Lu\ or a scalar i ?2 in the (3,2,1 /6) repre¬ 
sentation (the same as the quark doublet Q) that can couple to 
We note that dimension five operators involving 
these LQs can lead to proton decay; these have to be forbidden 
by a discrete symmetry, as shown in detail in 03. 

We focus on i ?2 because its quantum numbers are favorable 
from a DM perspective, as we will clarify later. Pair produc¬ 
tion of i ?2 via gluon fusion and quark-antiquark annihilation 
thus constitutes our main example. Its Lagrangian couplings 
are 

-A^^'4i7^eLi+h.c., (1) 

where i,j denote flavor indices. In the above. 



Expanding the SU(2) components yields 


- Ki daRiVaei - + h.c. ( 3 ) 

In addition to the LQ, we introduce a dark sector with a 
scalar S and a fermion %, with a Z 2 symmetry under which 
the dark sector is odd, whereas the SM and LQ sector is even. 
Thus, our new physics content is 

i?2 = (3,2,l/6)+ 

/ J-^-o ^-1- \ 

^ = ( 1 , 3 , 0 )_ = ( _^so j . (4) 

X = (1,1,0)-, (5) 

where we have also displayed the Z 2 quantum numbers as a 
subscript. Notice that there is almost no freedom in choos¬ 
ing the above particle content. The quantum charges of the 
dark sector are fixed by symmetries; either the scalar S or the 
fermion x has to be a triplet under SU(2 )l (having both sin¬ 
glets will not give enough charged leptons in the final state 
to match the CMS study, as we will see). We have chosen a 
scalar triplet for convenience; the discussion and results would 
be analogous for a model with a triplet fermion x and a sin¬ 
glet scalar S. The hypercharge of the dark sector is fixed to be 
zero to easily accommodate DM direct detection limits. 

The LQ decay into the dark sector can then be described by 
adding the following dimension-5 effective operators: 

- ^SQaR2 - ^Sl,xH , 

Ai A2 


where H is the iso-spin transformation of the Higgs dou- 
blet HSl. 

The CMS study does not discuss the relative number of 
e“''e“ and jjY events in the excess. We will 

assume them to be equal and introduce both a first generation 
LQ and a second generation LQ in order to obtain both e“''e“ 
and jjY events. In the following, however, we will describe 
our methods for first generation LQs, keeping in mind that 
second generation LQs can be similarly treated by a simple 
replacement of electrons by muons. 

Spectrum and Decays. Either the triplet component or 
the singlet x can be the lightest DM particle. The two pos¬ 
sibilities lead to distinct DM phenomenology. Here, we will 
mainly consider the case of a singlet x as the DM candidate. 

The two couplings of the LQ induce two decay 
modes:i ?2 —e/ and R 2 Sjx- The latter cou¬ 
pling in Eq. (§ induces the S ^ ex and S ^ vx decays of 
S. In components: —>■ xe^ , vx, so that de¬ 

cays of charged scalars give charged leptons and ]fjj in the 
final state. One loop electro-weak corrections induce a small 
mass splitting of « 200 MeV between the neutral and charged 
states in S. Combining all decays, one R 2 LQ can produce the 
following final states; 

1. A charged lepton and a jet. The free couplings of the 
Lagrangian allow us to set the i ?2 branching ratio to the 
level required to be compatible with CMS searches for 
first and second generation LQs, with final states eejj 
and etyjj, and similarly for muons. We will And that at 
the benchmark point of this study, there is no excess in 
these channels; thus, we will choose 

BR (^2 ^ ej) ~ 0%, 
consequently BR(i ?2 —>■ Sjx) ~ 100%. 

2. A jet and missing energy i^), with 

BR(i ?2 ^ S°jx ^ J^) ~ 29%. 

3. A charged lepton, a jet and missing energy, with 

BR(R2 ^ S^jx ^ e^Jpr) ~ 71%. 

Constraints on the masses of R 2 , S, x- (?) Constraints 
from jpT searches. LQ pairs decaying to jets and missing 
energy must be compatible with jfJr searches from ATLAS 
ETIl and CMS ll22ll . The most relevant exclusion limits from 
these searches are phrased in the squark-neutralino {rriq, rUj^o) 
mass plane (Eig. 10 of 11211 '). where the lightest neutralino is 
stable and escapes the detector. In our model, both 5° and 
X escape undetected. To compare to the experimental stud¬ 
ies, we map the neutralino mass m^o to ms + m^ and the 
squark mass mg to m^^. Taking into account that, in our 

case, BR(i ?2 —)■ jpr) ~ 29%, we conservatively estimate 
the bound 


( 6 ) 


ms + m^ > 300 GeV 


(7) 


3 


for LQs in the mass range 450 — 650 GeV. 

{ii) Constraints from CMS LQ searches: There is a mild 
evidence for ~ 550 — 650 GeV first generation scalar LQs, 
using 19.6 fb“^ of integrated luminosity, at 2.4 (t and 2.6a, 
in the eejj and ei/jj channels respectively i). A branching 
BR(i ?2 —>■ ej) ^ 15% agree well with the observed excess 
in the 550 — 650 GeV mass range. For second generation 
LQs, studies of the fifijj and iiiyjj final states has resulted 
in the exclusion of scalar LQs with masses below 1070 (785) 
GeV for = 1(0.5) Q. 

As stated before, we will assume a branching BR(i ?2 —i” 
ej) ^ 0% to switch off decays to this channel for LQs 
of mass around 500 GeV. We will also assume BR(ii 2 —> 
Hj) ^ 0% to be compatible with second generation LQ 
searches. 

{in) Constraints from the dilepton invariant mass mu 
distribution of the CMS excess, located mostly 

below mu ~ 80 GeV. The two leptons in our case come 
from the decay S £x- To get the excess in the required 
range the mass difference between S and x should be 
ms — m^ ~ 20 — 40 GeV . Spectra where ms — m^ > 40 
GeV are disfavored since the dilepton invariant mass distribu¬ 
tions would peak at a value of mu that is too large compared 
to the excess. On the other hand, for ms — m^ ^ 20 GeV, the 
leptons are too soft and do not survive the px cuts for leptons. 

Results. We take our background estimates from m. Op¬ 
posite sign opposite-flavor (OSOF) leptons from tt, which has 
the same rate of the same-flavor (OSSF) channel, are used to 
measure these backgrounds in the CMS study. Drell-Yan pro¬ 
duction, the main irreducible background, is estimated by a 
control region which does not overlap with the signal region. 

We follow the CMS counting experiment analysis in HI for 
the signal. The final state is required to have at least two lep¬ 
tons and at least two jets. The cuts employed are: 

Cut (i) Two OSSF leptons are required to be present, with 
px > 20 GeV in 1 77 1 < 1.4 which is defined as the central 
region in the CMS study. 

Cut (ii) At least two jets are required with px > 40 GeV in 
IpI < 3.0. 

We use FAST JET ll24ll to reconstruct jets. An event is se¬ 
lected if it contains two jets and satisfies pr > 150 GeV, or, if 
it contains three or more jets and satisfies pr > 100 GeV. In 
the dilepton invariant mass range 20 GeV < mu < 70 GeV, 
the total background estimate provided by the CMS study for 
central OSSF events is 730±40. The observed number was 
860, corresponding to an excess of ISO^t^g events provided by 
new physics, we hypothesize. In our model, this excess num¬ 
ber of events is produced by first and second generation LQs. 
We implement the model in FEYNRULES ll25]l and calcu¬ 
late the branching ratios and cross sections using MAD GRAPH 
5.11 ll26l . The events are then passed onto PYTHIA lIZTl 
for parton showering and hadronization followed by the mod¬ 
elling of detector effects by DELPHES 3.1.2 128]. We took 


exactly the same electron/muon/jets isolation criteria adopted 
in H]. Production cross sections were normalized by the NLO 
QCD rate using Prospino2.1 1291 . 

We performed a scan over the masses ms, frix)’ 

fixing the Lagrangian parameters Aj* = 10“^, hijAi = 10“^ 
GeV“^ and h{/ A 2 = 10 “^ GeV“^. In the above, we have 
assumed democratic values for the Yukawas corresponding 
to i = 1,2 and j = 1,2. We note that constraints aris¬ 
ing from meson decays and mixings are sensitive to prod¬ 
ucts of the Yukawas A®^; in particular, our choice of Yukawas 
is compatible with bounds given in 1^ and 1131 . The cut 
flow is displayed in Table |I] for some spectra that best fit the 
number of events in the central region of the CMS search 
with BR(i ?2 —>■ ej) ~ 0 %, ms + m^ > 300 GeV and 
ms — m^ < 40 GeV. 

We have assumed theoretical uncertainties following ED, 
which are the same as the estimates used by CMS in its LQ 
search ID. The uncertainties are given in Table 1 of the CMS 
paper ID. Variation of the renormalization/factorization scale 
between half and twice the LQ mass leads to a ~ 25% un¬ 
certainty in the production cross section. The PDF uncer¬ 
tainty is approximately 20% of the NLO cross section. This 
can easily lead to a ~ 20% variation in the number of sig¬ 
nal events, which will be interesting to pursue if the signal 
becomes stronger in the next run. 

In Fig. [T] we show the comparison between the CMS data 
points and the predicted distribution for our model at bench¬ 
mark point C. As we see, our model can fit the data very 
well. In fact, points A,B and C of Table [^present this feature. 
While the CMS study showed a triangular shape with a sharp 
edge, for benchmark point C around 21 events survive in the 
region mu > 100 GeV, which is well within the background 
uncertainty. 

In order to better quantify the agreement between these 
models predicitions and the data we used the x^ statistics 
defined below in Eq. which, as in the experimental study, 
compares a model prediction in the Central (\r]i\ < 1.4) and 
Forward (1.6 < < 2.4) regions 


Nbir, 

= mfo 

i—1 

+ + 

( 8 ) 

In this formula, and denote the i-th dilep¬ 

ton invariant mass bin for the signal, backgrounds(DY and 
FS) and data distributions, respectively, in the Central region. 
The superscipts I and O denotes the Inside the window re¬ 
gion 20 < mu < 70 GeV and the Outside region mu > 70 
GeV, respectively. The corresponding experimental uncertain¬ 
ties in each bin are given by cr/, a^ and a^. In order to 
better approach the uncertainties in the background bins we 
take into account two nuisance parameters 6 = ( 0 i, 02 ) for 
the DY background normalization, bi, and the FS background 


{psj +0-h\- dj)^ 
0.bO,-df)21 
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TABLE I. Summary of the effective cross-sections (fb) for some 
benchmark signal points that best fit the CMS signal at LHC8. In 
the fourth and fifth rows of each point we also show in parenthesis 
the final number of events predicted. The last row displays the num¬ 
ber of events in the forward region. Masses are in GeV. Lagrangian 
parameters are fixed at = 10“"^, hi/Ai — 10“^ GeV~^ and 
h'i/A 2 = 10“^ GeV~^. The benchmark point for the current study 
is C. The cross sections (in fb) predicted at LHC14 after imposing 
the cuts of the LHC8 analysis for the three best benchmark points is 
displayed in the last column. 



Selection 

Signal (fb) Signali 4 (fb) 


preselection 

149.1 

A : (450, 200,170) 

Cut (i) 

40.2 46.6 


Cut (u) 

10.5 


20 < m« < 70 

5.7(110) 


mu > 100 

1.8(37) 


Forward 

18 


B (500,200,170) 

preselection 
Cut (i) 

Cut (ii) 

20 < mu < 70 
mu > 100 
Forward 

105.7 

28.5 

11.6 

5.7(112) 

2.3(45) 

14 

69.6 


preselection 

126.1 


C : (500,160,140) 

Cut (i) 

19.0 

53.6 


Cut (ii) 

10.0 



20 < mu < 70 

5.8(114) 



mu > 100 

1.11(21) 



Forward 

8 



normalization 62 - 

All these regions are crucial to determine the pattern of de¬ 
cays and the spectra favored by the data. For example, those 
spectra where the ms — m^ are bigger than ~ 40 GeV are 
disfavored in this respect as their dilepton invariant mass dis¬ 
tributions peak outside the window. On the other hand, for 
ms — m^ < 20 GeV, the leptons are too soft and do not sur¬ 
vive the pt cuts for leptons. We show in Fig. ?? the dilepton 
invariant mass distribution for the point C. The CMS study 
compared their data against a R — parity conserving super- 
symmetric scenario where an sbottom decays to a SF lepton 
pair and missing energy. This type of decay presents a dilep¬ 
ton mass distribution with a triangular shape and a kinematic 
edge around m^ — m^o which fits the observed excesses well. 

The Forward region is an important control region since it is 
expected that heavy particles undergoing decays as in SUSY 
or LQ models (as suggested in this work) will produce central 
leptons. We checked that the number of events is small for 
the best LQ models, around 7 events for benchmark C and 
similarly for the others as shown in Table This is consistent 


with the CMS reported number of 6 ± 20 events in the forward 
region. Moreover, the expected fraction of signal events with 
3 or more leptons in those benchmark points is never beyond 
0.5%, again consistent with the absence of an excess in the 
trilepton channel. 

We show in Table the obtained with Eq. ® and the 
corresponding p-value calculated as p = 
where n = 60 degrees of freedom, for the LQ models, the 
SUSY model of the experimental study and the pure back¬ 
ground model. The numbers in parenthesis represent the one¬ 
sided probability in the tail of a Normal distribution. The 
CMS study assigns a 2.4cr significance for the SUSY model 
against 2.6cr from our estimate. The discrepancy might be 
due the estimation of the errors included in the fit which are 
not fully accessible to us. Nevertheless, we also found that 
the best fit is provided by a SUSY type signal and point B. 
The benchmark points A, B and C, by their turn, are the best 
fitting points in our model as can been seen in Table |I^ and 
very similar to the SUSY fit. 


Model 


p-value 

A 

35.0 

0.005(2.5(7) 

B 

34.3 

0.004(2.6(7) 

C 

35.0 

0.005(2.5(7) 

SUSY 

34.3 

0.004(2.6(7) 

Background 40.6 

0.03(1.8(7) 


TABLE 11. The statistics computed from Eq. for three benchmark 
LQ models, the SUSY model and the background model taken from the CMS 
analysis (T]. In the third column we show the p-value of each model and the 
corresponding statistical significance in parenthesis as described in the text. 

With the current data, none of these models can be excluded 
at 90% CL, for example, although some models can provide a 
better fit than others. The background model is currently the 
worst fit amongst all of them. Comparing the LQ and super- 
symmetry interpretations we see that the SUSY model is able 
to explain the data as well as our best point, B. The bench¬ 
mark point C, however, is the one that better fits the excess 
quoted in the CMS work in the signal region 20 < mu < 70 
GeV, 114 events against 130 of the CMS fit. Moreover, it pre¬ 
dicts the smallest number of events in the Forward region of 
all benchmark points, 8 events against 6 from CMS. 

In the last column of Table |I] we show the projected 
cross sections at the LHC 14 surviving the cuts of the LHC 8 
analysis for the three best benchmark points found. 

Comments on Dark Matter Phenomenology. DM stabil¬ 
ity is guaranteed by the discrete Z 2 symmetry. The options are 
to have either triplet or singlet DM that could be the fermion 
X or the scalar S as discussed thoroughly in ifTSll . Taking into 
account, direct, indirect and collider constraints the following 
possibilities emerge: 

1. A scalar singlet S coupled to the Higgs portal. After 
taking LUX, indirect detection and Higgs invisible de- 
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FIG. 2. The dilepton invariant mass distribution of the signal events 
(point C) is displayed. The vertical dashed lines indicate the invari¬ 
ant mass cut. Around 21 events survive in the region mee > 100 
GeV, which is well within the background uncertainty. 


cay widths into account the mass range ms > 100 
GeV and ms ~ 60 — 65 GeV is allowed. However, 
XENONIT is expected to masses up to a 1 TeV lf3^ . 

2. A fermionic singlet (x) with a significant pseudoscalar 
coupling to the Higgs portal results into a less con¬ 
strained spin-dependent scattering cross section ImI . 
With A being the scale of the portal interaction, and 
sin^ the pseudoscalar coupling, from Fig.7 of ll34l . it 
is clear that > 0(100) GeV is allowed by data for 
sin^ ^ ^ 0.7 with A ~ 1 — 5 TeV. The bounds do not 
depend much on whether x is Majorana or Dirac. 

Thus, a singlet fermionic DM candidate with mass 
m^ ~ 140 GeV that we have taken in our collider analysis of 
this paper is currently a viable option. 

Conclusions. In this paper, we have shown that the excess 
observed by CMS in the £'^i~jjpT search can be explained 
consistently within a class of LQ models. We particularly 
discussed a model which consists of scalar first generation 
and second generation LQs that decay dominantly to leptons, 
jets, and missing energy in the form of a stable DM candi¬ 
date. The confluence of proton decay constraints and DM di¬ 
rect/indirect detection results in a highly predictive model of 
LQs that can satisfy the CMS t^(.~ jjpr search. Our model 
provides a general proof of concept that a peak distribution in 
the required dilepton mass range can be obtained outside the 
purview of supersymmetry. 

Benchmark point C consists of first and second generation 
LQs with masses of 500 GeV, and dark sector particle masses 
of ms = 160 GeV and m^^ = 140 GeV. The LQs dominantly 


decay to (.jpr final states 71% branching), and subdomi- 
nantly into jfrr final states (~ 29% branching), while there is 
negligible branching into the canonical LQ final states of ij. 

While CMS fit a triangular shape in the opposite-sign- 
same-flavor dilepton invariant mass distribution after event 
selection and flavor subtraction, it is too premature to settle 
definitively on a kinematic edge. In our model, the number 
of signal events in the window 20 < mu < 70 GeV after 
event selection is 114 for the benchmark point. The dilepton 
mass distribution peaks in the window between 20 — 70 GeV 
with the required event count, while the number of events in 
the region mu > 100 GeV is within the background un¬ 
certainty. In a simple x^ fit. this point compares favorably 
with the CMS SUSY fit in 111. The model is consistent with 
the non-observation of signal in the forward region, and in 
the trilepton final states. We have also provided projections 
for the signal cross sections at LHC 14 in the best benchmark 
models. For the benchmark point C, for example, ~ 5360 
events are expected with 100 fb~^ luminosity. 

These results highlight the importance of refining the 
search strategy and that a potential LQ discovery is attainable 
in the next LHC run. 
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